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Toward Identification of Order Parameters in Skutterudites 
- a Wonderland of Strong Correlation Physics - 

Yoshio KuRAMOTO * and Annamaria Kiss ^ 
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Current status is described toward identifying unconventional order parameters in filled skutterudites 
with unique ordering phenomena. The order parameters in PrFe4Pi2 and PrRu4Pi2 are discussed in relation 
to associated crystalline electric field (CEF) states and angular form factors. By phenomenological Landau 
analysis, it is shown that a scalar order model explains most properties in both PrFe4Pi2 and PrRu4Pi2 with 
very different magnetic properties. In particular, the highly anisotropic susceptibility induced by uniaxial 
pressure in PrFe4Pi2 is explained in terms of two types of couplings. In the case of SmRu4Pi2, the main 
order parameter at low field is identified as magnetic octupoles. A microscopic mechanism is proposed 
how the dipole and octupole degrees of freedom mix under the point group T;, of skutterudites. 
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1. Introduction 

Filled skutterudite compounds provide an unprecedented 
framework where fundamental and long-standing problems in 
condensed-matter physics all show up within the same crys- 
] tal structure. In rare earth (R) skutterudites RT4X12 with T 
transition metals and X pnictogens, different combinations of 
constituent atoms lead to enormously rich variety of proper- 

■ ties, such as itinerant-localized dichotomy and hidden elec- 
tronic orders. At the present stage, theoretical work needs to 
combine a phenomenological Landau-type approach and mi- 
croscopic consideration for understanding the overall behav- 
ior. In this paper, we focus on identification of unconventional 
order parameters in PrFe4Pi2, PrRu4Pi2 and SmRu4Pi2. 

PrFe4Pi2 undergoes a second-order phase transition at To - 
6.5K, which accompanies a typical structure in the specific 

■ heat, a sharp peak in the magnetic susceptibility, and a steep 
increase of the resistivity just below the transition.^ A crystal- 
structure modulation with the wave vector Q = (1,0,0) was 
found below Tq by X-ray diffraction experiments,^ which 
is attributed to the existence of staggered local electronic 
states of the Pr ions. Early NMR^ and elastic measurements"^ 
were interpreted in terms of antiferro ordering of r3-type 
quadrupole moments. However, the persistent isotropy of 
the magnetic susceptibility in the ordered phase cannot be 
explained by r3 quadrupolar order. Furthermore, staggered 
dipoles are always parallel to the field direction both in neu- 
tron diffraction^'* and NMR.^ Recently, careful analysis of 
the NMR results have shown that the local symmetry at the 
Pr sites is preserved in the ordered phase. ^'^ Furthermore, the 
continuous field-angle dependence of the transition tempera- 
ture gives also evidence for the exclusion of r3 quadrupolar 
order.'*- 1" 

In a recent paper,"' we have proposed that the order param- 
eter in PrFe4Pi2 is a staggered electronic order which does not 
break the local Ti, symmetry around each Pr site. We call this 
order a scalar order since it has the Fig symmetry. It is found 
in ref. 10 that the scalar order scenario can explain naturally 
the isotropic magnetic susceptibility in the ordered phase, the 



field angle dependence of the transition temperature and mag- 
netization, and also the splitting pattern of the ^'P NMR spec- 
tra. In this paper, we proceed to CEF theoretical description 
of the scalar order, and compare PrFe4Pi2 and PrRu4Pi2. Fur- 
thermore, we explain not only the Neel-type anomalies of the 
magnetic susceptibility x in PrFe4Pi2 near the phase transi- 
tion, but also the huge anisotropy ofx induced under uniaxial 
pressure. 

As another prototype of mysterious orders, we take 
SmRu4Pi2, and analyze its CEF states. The conventional view 
point is that the CEF states are composed by linear combina- 
tion of Hund's rule ground state with J - 5/2. In this case, the 
highest rank of multipoles in this manifold is 2J - 5, which 
is too small to distinguish between the Oi, and Ti, symmetries. 
Namely, the Stevens operator of sixth rank O'^ - - O^, 
which makes the difference, ' ' has zero matrix elements with 
J - 5/2. On the other hand, recent experimental results sug- 
gest mixing of dipole and octupole degrees of freedom.'^' 
We analyze the wave functions in the CEF states taking higher 
order hybridization processes. It is found that the closeness of 
the 7 = 7/2 excited state above the 7 = 5/2 ground state tends 
to compensate the small ratio of hybridization over excitation 
energy of 4f^ intermediate states. 

2. CEF states and scalar orders in Pr skutterudites 

2.1 Relevant CEF states 

In previous work,'^ we ascribed the main source of CEF 
splittings to covalent hybridization effects between 4/ and 
ligand orbitals. The relevant point group Ti, causes mixing 
of two kinds of triplets r4 and Fs in the cubic case. The hy- 
bridized triplets are called F^'' with larger weight from r4, and 
F^^* with larger weight from Fs. For later use, we shall give ex- 
plicit form of these states in the case of 4/^ configuration with 
7 = L- 5= 5- l=4. 

|Fi>= 75/24(14) + 1-4)) V77T2|0>, (1) 
|F4;fl)= VT72(|4)-|-4)), (2) 
|F4;/7)= VT78|3)+ VtTsI-D, (3) 
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|F4;c)= Vl78|-3)+ VtTsU), (4) 
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|r5;fl)= Vl/2(|2)-|-2», (5) 

|r5;/7>= V778|3>- VT78I-1), (6) 

|r5;c)= V778I-3)- VTTSID, (7) 

in terms of eigenstates of J^. Similarly, the double Fa CEF 
states are given explicitly by 

\Ty, a) = 77/24(14) + |4)) - V57l2|0), (8) 

|r3;Z7)= VT72(|2> + |2». (9) 



The three representative Pr skutterudites, PrFe4Pi2, 
PrRu4Pi2, and PrOs4Sbi2 all have the singlet CEF ground 
state. However, the first excited level is different from each 
other. Namely, PrFe4Pi2 has the low-lying F^' ' triplet with 
strong van-Vleck susceptibility, and possibly the F23 dou- 
blet, which goes over to F3 doublet in the cubic symmetry. 
We suspect that these six levels are almost degenerate in the 
high-temperature phase of PrFe4Pi2. As temperature becomes 
lower than To, one of the Pr sublattices has a singlet CEF 
ground state, while the other Pr appears to take the doublet. 
This conjecture comes from neutron scattering of PrFe4Pi2 
where at least two inelastic transitions are visible in the or- 
dered phase. 

On the other hand, PrOs4Sbi2 has the low-lying F^^' triplet 
which goes over to F5 in the cubic point-group symme- 
try. Hence the quadrupolar (van-Vleck) susceptibility is large 
in PrOs4Sbi2. Finally, the triplet in PrRu4Pi2 in the high- 
temperature phase appears to be a strong mixture of F4 and F5 . 
In the ordered phase, one of the Pr sublattices has the singlet 
CEF ground state, while the other has the crossing of singlet 

(7.) 

and F4 triplet levels with decreasing temperature. Further de- 
crease of temperature brings about the point group lower than 
Th as observed in the splitting of the triplet of the order of 1 

2.2 Angular form factors associated with the scalar order 

Multipolar interactions of rank four (hexadecapole) or rank 
six (hexacontatetrapole) have a chance to bring an electronic 
order which keeps the original T/, symmetry around each Pr 
site, but lead to A and B sublattices with different CEF ground 
states. The shape of the Fermi surface with good nesting prop- 
erty should be responsible for the staggered AB sublattice 
structure.'^' Since the scalar order accompanies a slight lat- 
tice distortion, it can be probed by X-ray diffraction.^ More 
detailed information should be obtained if 4/ form factors are 
probed by azimuthal scan in resonant X-ray scattering using 
the electric quadrupole (£2) channel.'"^ 

Since the fourth-rank tensor relevant to E2 scattering is 
very complicated, we visualize the scalar order by deriving 
the simplest form factor that corresponds to a weighted av- 
erage of the electron charge density. Namely, we utilize the 
integer (7 = 4) value of the Pr^^ configuration, and introduce 
a fictitious "wave function": 

ilJra{^)^{D.\T,a), (10) 

where Q. represents the solid angle specified by {0, 0) such 
that c/Q = sm0d6d(p. Then iffrai^) can be derived in terms 
of spherical harmonics Y4„{Q.) with use of eqs.(l) to (9). The 
angular form factor prjiO.) associated with a CEF level F is 
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Fig. 1. Angular form factor of Fi singlet CEF state. 

defined by 

pry(n) = ^w„|iAr«(S^)l', (H) 

a 

where Wa is the weight factor of the component a in the CEF 
states. In the case of singlet Fi, we have Wq. = 1, while in 
the case of doublet or triplet, we have Wa - 1/2 or Wa - 
1/3, respectively. Since IF, a) is not a single-particle state, and 
since spins are also involved, the form factor prj{i^) is not a 
charge density itself. One may nevertheless gain good insight 
into spatial pattern of the scalar order by prjiQ.). 

Figure 1 illustrates the angular form factor associated with 
the CEF singlet state Fi . The distance of a point on the surface 
from the origin represents pr y(f2) for each solid angle. Clearly 
the cubic symmetry is preserved in the form factor, which is 
understood as a superposition of a constant, a hexadecapole 

+ y'^ + - 3/5 with r = (x, y, z) being a unit vector, and a 
hexacontatetrapole (x^ - y^){y^ - z^)(z^ - x^). In the case of 
the multiplet, on the other hand, each wave function breaks 
the cubic symmetry. However, summation over equally dis- 
tributed degenerate states recovers the pattern consistent with 
the cubic symmetry. 

Figure 2 shows the form factor of each F3 state (a) and 
(b), and the average of the two in (c). Figure 3 shows the 
scalar form factor associated with the F5 triplet (a), and the 
r4 triplet (b). It is seen that F5 has a pattern dominated by a 
hexadecapole, while the patterns of Fj, and F4 are dominated 
by a hexacontatetrapole. In the case of PrFe4Pi2, we expect 
the form factor in Fig. 1 is realized in the A sublattice, while 
the form factor in Fig. 2(c) in the B sublattice. This is the mi- 
croscopic image of the scalar order in PrFe4Pi2. On the other 
hand, in the case of PrRu4Pi2, we expect the form factor in 
Fig.l is realized in the A sublattice, while the form factor in 
Fig. 3(a) in the B sublattice. 

Figure 4 illustrates the staggered arrangement of angular 
form factors on the bcc lattice formed by rare-earth ions in 
skutterudites. The spatial symmetry in the ordered phase re- 
mains cubic, but the unit cell is doubled. As a result, the su- 
perlattice has a simple cubic structure. 

As temperature approaches to zero, one of the sublattices 
has to release the nonzero entropy associated with the dou- 
blet or triplet degenerate level. In PrRu4Pi2, the triplet seems 
to split into a singlet and a doublet by about 1 K.'^ On the 
other hand, there is no information about the CEF states in 
the ordered phase of PrFe4Pi2. If one of the sublattice is the 
doublet, either the lattice distortion or the quadrupolar Kondo 
effect should break the degeneracy. Such low-temperature be- 
havior deserves further experimental study. 
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Fig. 2. Angular form factor of doublet states: (a) IPs ; a) state; (b) IFj ; h) 
state; (c) average over the both states. 





Fig. 3. Angular form factor of triplet states averaged over three compo- 
nents: (a) Fs state; (b) F4 state. 



3. Landau expansion of the free energy 

We use the phenomenological description of coupling be- 
tween the scalar order parameter and other degrees of freedom 
such as magnetization, quadrupole moment, and lattice strain. 
In the Landau theory, one expands the free energy in terms of 
the set of electronic order parameters "P, , which are taken to 
be real. Up to fourth-order, we write 



1 



_f...vi<2m2 



(12) 



where we have introduced the quantity a, - UjiT - T,). The 
constants a,, fe, are positive, while gij and c,j can have either 
sign, r, is a hypothetical transition temperature without cou- 
pling to other order parameters. The actual transition occurs 
at To corresponding to the scalar order, which gives the largest 
of all r,. For other component *P,, we neglect the correspond- 
ing bi in most cases. These parameterizations have a merit 
that each coefficient is regarded as a constant as long as the 
temperature is close to Tq, and the external perturbations are 
small. 

As explicit constituents of ^F,, we include the scalar or- 
der parameter i/tq, the homogeneous magnetization M, the 
r3-type homogeneous quadrupoles. Furthermore, we also in- 
clude the lattice strain components Sxx, Syv, Szz, which have a 
bilinear coupling with quadrupole moments Qij. The second- 




Fig. 4. Illustration of the sublattice form factors of the scalar order: (a) Fi- 
Fj staggered order; (b) F1-F3 staggered order. 



order couplings C/e with / - Q, ij/Q can be neglected because 
the background elastic constant cf.' is large enough. We de- 
fine the quadrupole moments by 

Qu = 0° = (1/ V6)(272 - j] - 4), (13) 
Q,.^Ol^{\l^){Jl-Jl), (14) 

and introduce the notation - Ql + Ql- Similar notations 
e„, St, are used for the strain components. Then the magneto- 
elastic couplings are described by 



(15) 



wherefij - (1/ V3)(£x.f + eyv+£zc)- In addition, the free energy 
in eq.(12) includes the coupling 

gMQ ^ (2M2 -Ml- Ml) Qu + ^ [Mj - Ml) e, 

which we call the tensor coupling. 
4. Magnetic susceptibility 

4.1 Coupling of moments with the order parameter 

The magnetic susceptibility is obtained from the formula: 
- d^T IdM^. We first consider the case without uniaxial 
pressure. For T > Tq, we set i/'g 
Weiss law: 



0, and obtain the Curie- 



X+ = aM(T - Tf), 



(16) 



where Tf is the Weiss temperature, and the superscript + in- 
dicates T > To. In the ordered phase close to the transition 
temperature, we obtain 



xJ = aM(T - Tf) + c^M^Q + ^gMsSs- 



(17) 



Since the system is far from the volume collapse, the last 
term with gMs is almost a constant, and can be incorporated 
into renormalization of Tf. Using the equilibrium condition 
d'F/diffQ = 0, we eliminate i/Tq and obtain the inverse mag- 
netic susceptibility as 
-1 



X- 



auiT - Tf) - c^Ma^iT - To)/b^. 



(18) 



There is a peak in the magnetic susceptibility at T - To if 
dxZ^ /dTlfa < is satisfied. The temperature derivative is cal- 
culated as 



dT 



— aM - 



(19) 



which gives the condition a^c^M > anbip for occurrence of 
the peak \nx{T) at T = Tq. The peak means that the growth of 
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the order parameter gives a negative feedback to the magnetic 
fluctuation through the repulsive coupUng c^m- 

Expressions (16) and (18) for the magnetic susceptibility 
are used to fit the measured result. We set Tp = 3.5K for 
the Weiss temperature,' and value for um is obtained from 
the fit in the paramagnetic phase. The value for the combina- 
tion a^CfMlbf is obtained from the fit of the susceptibility in 
the ordered phase. The parameters b^, and c^m are further 
constrained by the experimental temperature-magnetic field 
phase boundary. The fitting with account of these constraints 
is shown in the upper panel of Fig. 5. 

In the case of PrRu4Pi2, the anomaly ofx at the scalar tran- 
sition is very small. We interpret the smallness in terms of the 
foUowings: (i) the system is far from the ferromagnetic insta- 
bility characterized by Tp in eq. (16), and (ii) the couphng 
constant c^m is small. In fact the elastic anomaly at the tran- 
sition temperature Tq ~ 70K is less pronounced as compared 
withPrFe4Pi2.2° 

4.2 Uniaxial pressure effect 

Intriguing experimental results are obtained for the mag- 
netic susceptibility in the presence of uniaxial stress.-^'' ■^•^ 
Namely, the magnetic susceptibility shows large enhancement 
for the uniaxial pressure applied parallel to the magnetic field 
direction iH\\cr), while it shows only slight decrease when 
the pressure is applied perpendicular to the field direction 
(H ± cr). 

Now we discuss the properties of the magnetic susceptibil- 
ity around the transition temperature Tq. The direction of the 
uniaxial stress is taken as cr||(001), and we consider two dif- 
ferent directions of the magnetic field, namely /f||(001) and 
H\\{100). For small values of the uniaxial stress, it is enough 
to consider only the linear term in cr. Thus, we obtain the sus- 
ceptibilities Afn for //||(001) andx± for //||(100) as 

-1 2 4 

=aM + c^M^I^Q + IgMsSs + -^uqQu, (20) 

2 

X'^ = ttM + C^iM^Q^ + IgUeSs -pgMQQu- (21) 

V6 

Among the terms appearing in the expression of the magnetic 
susceptibility, gMe^s and gMoQu contain the uniaxial stress. 
We find in eqs.(20) and (21) that the former term (scalar) is 
isotropic, while the latter term (tensor) is anisotropic and has 
different sign for the two magnetic field directions. There- 
fore, the cancellation of the scalar and tensor terms can oc- 
cur for the case H ± cr, which reproduces the experimen- 
tal situation. Now we set t/zQ, e., and Qu from the equilib- 
rium conditions d7^/di//Q = 0, dT jdss = (Ts = VT73o" and 
&Fjdeu = (Tu = V2/3cr. Keeping only the leading term as 
in calculation without uniaxial stress, we obtain the magnetic 
susceptibihties^ll and;ifj_ in the ordered phase as 



^11. 
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where x- is given by eq.(18). Here the relevant elastic con- 
stants are written as C3 = cf^ - cf^, Co = cf^ + 2Cf^, and 
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Fig. 5. Top: Magnetic susceptibility around the transition temperature Tq. 
The parameter values are = lO^'Pa], Tf = 3.5K, om = 9.6 ■ 
103[Pa-/jg2.K-l], = 1.95 ■ lO'^Pa-R-'], c>m = 1.37 ■ lO'^Pa -pg^]. 
Boxes represent the measured result taken from ref. 1. Bottom: Boxes rep- 
resent the measured result taken from ref. 21. 



The susceptibilities in the paramagnetic phase can be obtained 
by taking = in the eqs.(22) and (23). Experimentally 
magnetic susceptibiUty at the transition temperature depends 
linearly on the uniaxial pressure for small values of cr. With 
a proper choice for the values of gMQ and gMe, we can fit 
the measured susceptibilities ;k'||(7'c) and;^'j_(7'e) at the transi- 
tion temperature. The result can be seen in the lower panel 
of Fig. 5. We consider that the almost constant behavior of 
the magnetic susceptibility ;^fx(7'c) at T = Tc is accidental. 
Namely, different anisotropic behaviors are also possible de- 
pending on the parameters gMs and gMQ- 

5. Mixing of different angular momenta in CEF states 

The CEF splitting is caused not only by aspherical 
charge distribution around each rare-earth site, but also by 
anisotropic hybridization processes.'^ The hybridization is 
taken in the form: 



(25) 



Tver 



(24) 



where T is an irreducible representation in T/,, and v is an 
element therein. If there are different electronic states with 
the same representation, we distinguish them in terms of the 
index a such as r(a). 

In the standard theory for rare-earth, the CEF states are ob- 
tained by diagonalizing the Hund's rule ground states with 
given / under the CEF potential. The mixing of excited states 
with different J is neglected because the spin-orbit splitting 
is much larger than the typical CEF splitting. However, the 
difference of the energy splittings alone does not justify the 
neglect of higher J because the degree of mixing is not only 
determined by the Coulombic CEF potential, but also by the 
covalent hybridization with much larger energy scale. In order 
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4/'^ + ligandhole 




J =5/2 



Fig. 6. Second-order perturbation processes to mix J = 7/2 states in the 
CEF ground state. The spin-orbit splitting Ago is about 0.12 eV in Sm^*. 



to clarify the situation, we use the BrUlouin-Wigner perturba- 
tion theory which gives the formally exact series of perturbed 
wave function in terms of unperturbed on O as follows: 



»P = (D-H 



Q 



E-Ho 



Q 



E-Ha 



//hyb (26) 



where Q is the projection operator to make states orthogo- 
nal to <I>, and E is the exact energy of *P. The unperturbed 
Hamiltonian Hq describes the decoupled f-electron and lig- 
and states. Then the intermediate state i/hyb'l' is dominated by 
4yn-i (4y"+i ) states plus an extra ligand electron (hole). In 
the following we take the specific case where the trivalent Sm 
with 4/^ configuration is the ground state, which has domi- 
nant intermediate states with 4/*". The weight of the (9(//hyb) 
term in ^' is small because the hybridization is less than 1 eV, 
and the excitation energy to Sm^^ is typically a few eV. 

On the other hand, some the 0{H^^^ term has a larger 
weight in *P than the first-order term. This is because the sec- 
ond intermediate states are dominated by the Sm-'^ configu- 
ration which has the same 4/ electron number as the ground 
state but with different values of J. Then the excitation en- 
ergy in £ - //() is of the order of spin-orbit splitting, which 
is especially small in the case of 4/^ with merely 0.12 eV. 
If we put ffhyb ~ 0.3eV and f/eff ~ E - Hq ~ 2eV, we can 
roughly estimate as H^^^HUesh^o) ~ 0.4. Hence the ratio 
of hybridization over the excitation energy exceeds unity, and 
the weight of mixed states with 7 = 7/2 or higher is not neg- 
Ugible. Figure 6 illustrates the perturbation processes. 

The ground CEF states have wave functions with / = 
5/2: 

■ ~ (27) 

(28) 



\a±) = V5/6I ± 5/2) + Vl/61 + 3/2), 
l^'i) = I ± l/2>, 



where a and b specify the orbital quantum number, and ± the 
Kramers partners. In the case of 7 = 7/2, the wave functions 
with Fgv symmetry are given by 



|a±) = VlMI ± 5/2) + V374I + 3/2), 



\I3±) = ^JVT2\ ± 7/2) - V57T2I + 1/2), 



(29) 
(30) 



where the orbitals are now specified by a,/3. Any linear com- 
bination of a 7 = 5/2 state in Fgv and that with 7 = 7/2 sat- 
isfies the symmetry requirement for the group. Hence the 
eigen function of the actual CEF potential is determined by 
the mixing term connecting 7 = 5/2 and 7/2 manifolds. Ac- 
cording to our estimate in the preceding paragraph, the mag- 



nitude of the effective mixing potential 

Q 



Veif = -f^hyb 



E-Ho 



Hhyb 



(31) 



should be larger than, or at least of the same order of mag- 
nitude as the spin-orbit splitting Aso ~ 0.12 eV. If the Sm^^ 
state is closer to mixed valence, the excitation energy to 4/^ 
(Sm^"^) should also be small. As a result, the CEF ground state 
is expected to have considerable weight of 7 = 7/2 and higher 
angular momenta. It should be interesting to estimate the mix- 
ing more quantitatively. 



Bilinear coupling 
SmRu4Pi2 



of dipoles and octupoles in 



The most important consequence of the mixing of 7 = 7/2 
and even larger angular momentum in the Fg? CEF state is 
that the matrix element of the sixth-rank tensor 01 = Ol - 

6 6 

(9g becomes nonzero. Then the symmetry Tij lower than (9/, 
makes mixing of two triplet representations r4 and Fs, which 
has been recognized important in Pr skutterudites. In the T/, 
symmetry, there is a single triplet representation called F4, 
and different states are specified by a superscript such as F^ ' 

(2) 

or . Physically speaking, most dipoles and octupoles mix 
in the T/, group. The only pure octupole is the pseudo-scalar 
Fi„ which transforms as a third-rank tensor xyz. In Oh, the 
odd representation F5„, which transforms as x(y^ - z^) and its 
cyclic partners, is also a pure octupole. 

It has been suggested by Yoshizawa et al. that the ordered 
phase in SmRu4Pi2 has the dominant F5„ octupole compo- 
nent, which mixes with the dipole component r4„ in the nota- 
tion of Oh- Yoshizawa' s idea is most simply illustrated by the 
following form of the Landau free energy: 



1 



1 



T =-a,{T-T^)% + -as{T■ 



Ts)% + VY4*I'5, (32) 



with neglect of higher order terms. The last term with a real 
coupling constant v represents the mixing between dipoles 
and octupoles represented by *I'4 and *P5, respectively. With- 
out the coupling term, each order would have set in at tem- 
perature 7^4 or T^. The actual transition temperature is given 
by 

..2 ll/2 



7^5)' + 



(33) 



while the partner temperature T^. with the negative sign for 
the square root in eq.(33) is not a true transition, since the 
order parameter already grows below 7^+. If the coupUng v 
is small, however, Tc- may appear as a crossover with some 
structure in physical quantities. 

This scenario explains most naturally the appearance of the 
internal magnetic field,^^"^^ ferro-quadrupole moment,'^ and 
the large elastic anomaly .^^-^^ These features are analogous 
to those in the phase IV of Cei-xLaxB^ with x ~ 0.7, where 
the antiferro octupole order has been proposed,^^'^^ and con- 
firmed by various measurements. ^''"^^ The cubic symmetry in 
Cei_;cLa;cB6 allows pure octupoles which should have zero 
internal field at the Ce nucleus site, but with finite off-center 
field which has recently be detected by neutron scattering.-'^ 

The domain structure of SmRu4Pi2 which is consistent 
with the NQR result^^ is the same as that in Cei_;cLa;(:B6; there 
appear ferroquadupole domains with principal axis along 
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Fig. 7. Preferred quadrupole domain for magnetic field in the xy-plane with 
<9 < n/2 (Left). In magnetization and elastic constants in SmRu4Pi2, a 
domain switching to [1 1 1] for n/2 <6 <n should lead to a cusp structure 
on top of the two-fold pattern shown by the dotted line (Right). 



[Ill] and other three equivalent directions. If a magnetic field 
favors a quadrupole whose principal longer axis is the clos- 
est to the field direction, a sudden switching from one do- 
main to another should take place as magnetic field is rotated. 
This switching should appear as a cusp structure in physical 
quantities such as magnetization and elastic constant. In fact, 
Yoshizawa et al. have recently observed an intriguing pattern 
in the elastic constant of SmRu4Pi2.^'* Figure 7 illustrates our 
interpretation of the pattern. 

Experimentally, the lower transition becomes more and 
more visible as the applied magnetic field becomes stronger. 
In the above scenario, it seems difficult to reproduce such a 
behavior as long as the couphng v is independent of mag- 
netic field. Another possibility is that the second transition is 
also a real phase transition above a critical magnetic field, but 
becomes a crossover below the critical field. Then the pres- 
ence of critical point at finite temperature and field requires 
another model which is certainly more complicated than the 
one given by eq.(32). Within the Landau phenomenology, we 
have been unable to find a model with the desired property; 
the second transition, if real, always starts from zero tempera- 
ture instead of a critical point at finite temperature. It remains 
a challenge for theory to construct a microscopic model which 
goes beyond the scope of the Landau phenomenology. From 
experimental side, it is desirable to decide whether the second 
transition is a sharp crossover, or a real transition. 
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